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ABSTRACT 

The protostellar mass function (PMF) is the Present-Day Mass Function of the protostars in a 
region of star formation. It is determined by the initial mass function weighted by the accretion time. 
The PMF thus depends on the accretion history of protostars and in principle provides a powerful 
tool for observationally distinguishing different protostellar accretion models. We consider three basic 
models here: the Isothermal Sphere model (Shu 1977), the Turbulent Core model (McKee & Tan 
2003), and an approximate representation of the Competitive Accretion model (Bonnell et al. 1997, 
2001a). We also consider modified versions of these accretion models, in which the accretion rate 
tapers off linearly in time. Finally, we allow for an overall acceleration in the rate of star formation. 
At present, it is not possible to directly determine the PMF since protostellar masses are not currently 
measurable. We carry out an approximate comparison of predicted PMFs with observation by using 
the theory to infer the conditions in the ambient medium in several star-forming regions. Tapered and 
accelerating models generally agree better with observed star-formation times than models without 
tapering or acceleration, but uncertainties in the accretion models and in the observations do not 
allow one to rule out any of the proposed models at present. The PMF is essential for the calculation 
of the Protostellar Luminosity Function, however, and this enables stronger conclusions to be drawn 
(Offner & McKee 2010). 

Subject headings: stars: formation stars: luminosity function, mass function 
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1. INTRODUCTION 

The Initial Mass Function (IMF) is of central impor- 
tance in star formation since the properties of a star and 
its effects on the surrounding medium are determined pri- 
marily by its initial mass. One of the main approaches 
for inferring the IMF is to apply evolutionary models to 
observations of the Present-Day Mass Function (PDMF) 
of a group of stars. The IMF is created during the process 
of star formation, when the mass of each protostar grows 
by accretion until it reaches its final value. Observations 
of a region of star formation can, in principle, allow one 
to infer the PDMF of the protostars there; we refer to 
this as the Protostellar Mass Function (PMF). The PMF 
depends on both the IMF, which determines the relative 
number of stars when the star formation is complete, and 
on the process of star formation, which determines how 
long each protostar spends at a given mass. Because of 
this latter property, the PMF is potentially a powerful 
tool for inferring the nature of the star formation process. 
For example, inside-out collapse of an isothermal sphere 
(Shu 1977) leads to protostellar lifetimes that are propor- 
tional to the mass of the final star, whereas models based 
on competitive accretion (Zinnecker 1982; Bonnell et al. 
1997) have protostellar lifetimes that are independent of 
the mass of the final star; the inside-out collapse of a 
turbulent core (McKee & Tan 2002, 2003) has an inter- 
mediate behavior. As we shall see below, it follows that 
the PMFs predicted by these models are quite different. 

There are of course other approaches for observation- 
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ally distinguishing among the different theories of star 
formation. One is to study the relation between the 
mass distribution of density concentrations in molecu- 
lar clouds (the Core Mass Function, or CMF) and the 
stellar IMF, which are observed to be similar (McKee & 
Ostriker 2007). This similarity is the basis for recent the- 
ories of the IMF, which are predicated upon the assump- 
tion that stellar masses are determined by the produc- 
tion of gravitationally bound cores in turbulent molecu- 
lar clouds (Padoan & Nordlund 2002; Padoan et al. 2007; 
Hennebelle & Chabrier 2008, 2009). Such a direct con- 
nection between the CMF and IMF would appear to be 
a natural prediction of theories based on the inside-out 
collapse of gravitationally bound cores (Shu 1977; Mc- 
Kee & Tan 2002), but inconsistent with the theory of 
competitive accretion (e.g., Bonnell et al. 1997). Indeed, 
in inside-out collapse theories, the PMF and CMF are 
closely related, since protostars are forming in a signifi- 
cant fraction of cores. Clark et al. (2007) have pointed 
out that the mass distribution of cores depends on their 
lifetime, just as we shall see below that the PMF de- 
pends on the accretion time, and that this dependence 
would make the slope of the IMF significantly steeper 
than that of the CMF. However, in the Turbulent Core 
model (McKee & Tan 2003) the star formation time de- 
pends only weakly on the mass of the core, and (Hen- 
nebelle & Chabrier 2009) have argued that as a result the 
slopes of the CMF and IMF would agree within the ob- 
servational errors. It is clear that additional approaches 
for testing theories of star formation are needed. 

The principal difficulty with the PMF method devel- 
oped here is that at present, it is difficult to infer the 
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masses of individual protostars, both because the evolu- 
tionary tracks of protostars are uncertain (e.g., Chabrier 
et al. 2007) and because of the effects of accretion on the 
spectrum of the protostar, which are difficult to quantify. 
These difficulties will presumably be overcome in the fu- 
ture. Currently, the best way to determine the PMF 
appears to be through observations of the Protostellar 
Luminosity Function (PLF), which will be discussed in 
Offner & McKee (2010; hereafter Paper II). 

In this work, we define a protostar as an embedded 
source that is still experiencing signficant accretion and 
thus is more than a few percent from its initial stellar 
mass. Observationally, young stellar objects are charac- 
terized on the basis of the slopes of their spectral energy 
distributions into four classes, O-III (Adams et al. 1987; 
Andre & Montmerle 1994). However, since geometric ef- 
fects due to disk and outflow orientation influence the 
radiation reprocessing of the embedded source, it is dif- 
ficult to directly map classes to physical stages. For the 
latter, we use the classification proposed by Crapsi et al. 
(2008): Stage objects have protostellar masses that are 
less than or equal to their envelope, while Stage I objects 
have protostellar masses larger than the envelope. Once 
the envelope mass falls below 0.1 Mq, the object is con- 
sidered to have completed its main accretion phase and 
entered the Stage II phase. Although these stages do not 
directly correlate with class definitions, in practice. Class 
and Class I sources approximately correspond to Stage 
and Stage I (see §6 for additional discussion). 

In our analysis, we shall make two main assumptions: 
First, we shall generally adopt a Chabrier (2005) IMF 
with an imposed upper mass cutoff at m^,. This IMF 
has a log-normal form below \Mq and a power-law form 
above. For star formation regions large enough to fully 
sample the IMF, m„ ~ 150Mq (Figer 2005). However, 
as we show below, regions of low-mass star formation of- 
ten have maximum stellar masses of only a few Mq , even 
though in some cases one would expect stars more mas- 
sive than that according to the Chabrier (2005) IMF. It 
does not appear that this is a selection effect, since the 
data we compare with (Evans et al. 2009) represent a 
complete survey of nearby star-forming regions. If the 
deficit of high-mass stars is not a selection effect or a 
rare statistical fluctuation, then it must be the result 
of a physical inability to form high-mass stars in some 
regions — e.g., Krumholz & McKee (2008). In any case, 
our assumption is that the IMF above \Mq in these re- 
gions can be approximated by a truncated power-law im 
which the upper limit on the final masses of the proto- 
stars is set by the inferred upper limit on the masses of 
the more numerous newly formed stars (primarily Class 
II sources) in the sample. In the applications to observa- 
tions below, we shall generally adopt an upper mass limit 
of ZMq. Second, we shall assume that the accretion rate 
onto the protostar can be expressed as a simple function 
of the instantaneous protostellar mass, the final proto- 
stellar mass (i.e., the initial stellar mass), and the time. 
In particular, we ignore complications associated with an 
initial Larson-Penston accretion phase, when the accre- 
tion rate can be much larger than the value expected on 
the basis of dimensional analysis (e.g., McKee & Ostriker 
2007; however, in one of the most complete simulations 
of the formation of a star to date, Machida et al. (2009) 
found only a small enhancement in the accretion rate at 



early times); Furthermore, we average over any tempo- 
ral fluctuations in the accretion rate, such as occur in 
FU Ori outbursts (Hartmann & Kenyon 1996). In prin- 
ciple, the accretion rate can depend on the location of 
the protostar in its natal cluster; in our treatment, any 
such dependence is encoded in the dependence on the 
final mass of the protostar. 

The protostellar mass and luminosity functions were 
first considered by Fletcher & Stabler (1994a, b). In their 
two companion papers, they derived the time-dependent 
mass and luminosity functions for young embedded stel- 
lar clusters, including the luminosity contributions from 
protostars through main-sequence stars. Our work dif- 
fers from theirs in several important respects: (1) We 
consider a variety of different theories for star formation, 
allowing for the accretion rate to depend on the proto- 
stellar mass, m, and the final stellar mass, my, whereas 
they assumed that the accretion rate for all the stars was 
constant. Non-constant accretion allows us to test dif- 
ferent theories of star formation. In particular, as noted 
by Shu, Adams, & Lizano (1987), the constant accretion- 
rate model, which was developed for low-mass star for- 
mation, is unlikely to apply to high-mass star forma- 
tion. (2) As we mention above, observations in the solar 
neighborhood suggest an upper cutoff in the mass distri- 
bution at a few solar masses; we explicitly allow for this 
in our adopted IMF, whereas Fletcher & Stabler did not. 
(3) They treated the protostellar, pre-main sequence and 
main sequence stages, whereas we consider only the pro- 
tostellar case. (4) They focused on the case in which the 
star formation rate is constant for a given time period, 
whereas we consider both the steady-state case and the 
case of accelerating star formation (Palla & Stabler 1999, 
2000). (5) Finally, we note a difference in approach: they 
determined the probability that a star of a given mass 
would cease accreting and become a pre-main sequence 
star, whereas we label each protostar by its instanta- 
neous mass and its final mass. As we shall see below, 
this makes it straightforward for us to consider the case 
of tapered accretion, in which the accretion rate declines 
prior to the end of the protostellar stage. 

We begin with the derivation of the PMF in terms of 
the IMF and the accretion history of the protostars in §2. 
We consider four accretion histories in §3: the collapse 
of an isothermal sphere (Shu 1977), the Turbulent Core 
model (McKee & Tan 2002, 2003), a blend of these two 
models, and the Competitive Accretion model (Bonnell 
et al. 1997). We also consider the effects of tapering the 
accretion rate as the mass approaches the final mass, as 
suggested by Myers et al. (1998). In §4, we evaluate the 
PMF for these accretion histories, both analytically for 
a Salpeter (1955) IMF and numerically for a Chabrier 
(2005) IMF. PaUa & Stabler (2000) have suggested that 
nearby low-mass star-forming regions have star forma- 
tion rates that are accelerating in time, and we show how 
this affects the PMF in §5. We make a brief comparison 
with observation in §6 (a more extensive comparison will 
be given in Paper II) and summarize our conclusions in 
§7. 

2. STEADY-STATE PROTOSTELLAR MASS FUNCTION 

Consider a region in which stars are forming. We shall 
refer to this group of stars as a cluster, although we make 
no assumption as to whether the group of stars is grav- 
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itationally bound. In this section, wc assume that star 
formation is in a steady state, so that there are both 
stars that have reached their final mass and protostars. 
For clarity, we list our mathematical symbols and their 
definitions in Table 1. Each protostar begins with a neg- 
ligible initial mass and accretes until it reaches a final 
mass nif. The cluster IMF describes the distribution of 
stars in the cluster with respect to their final mass: in 
a range of final masses druf, the rate at which stars are 
forming is 

cLA/"* = A/'*^(TO/)(ilnm/, (1) 

where A/"* is the total star formation rate (i.e., the num- 
ber of stars forming per unit time). In the cluster, stel- 
lar masses extend from a lower limit me (theoretically 
expected to be O.OO4M0^Low & Lynden-Bell 1976 
as updated by Whitworth et al. 2007) to an upper limit 
m„, which we infer from observations of the cluster. Note 
that the IMF is normalized to unity, 

V'(m/)dlnm/ 



1. 



(2) 



As discussed in §1, one of our assumptions is that, in 
a steady state, the accretion history of a protostar is de- 
termined by two parameters, its mass and its final mass. 
The distribution of protostars is therefore completely de- 
scribed by the bivariate PMF, tpp2{m, w/), such that the 
number of protostars in the mass range dm with final 
masses in dmj is 

d^Afp = N'ptpp2{Tn,mf)dhim dlnmf, (3) 
where Afp is the total number of protostars in the cluster. 
The normalization of the bivariate PMF can be expressed 
in two forms: 

dlnmf / d\um ipp2{m,mf) = 1, (4) 



/ rflnm / dlnmf ipp2{m,mf) = l, 

Jo Jruf^e 



(5) 



where the ranges of integration explicitly impose the con- 
dition that mf > m and where 

mf^i=max{m,mi) (6) 
is the lower bound on the integration over m/ in equation 
(5). The region of integration in the m — mf plane is 
shown in Figure 1. 

The PMF is the present day mass function of the pro- 
tostars in the cluster. 



il)p{m) = i ipp2{m,mf)dlnmf; 

Jrrif^e 



(7) 



it is normalized to \inity also, as ensured by equation 
(5). The PMF is observable in principle, although that 
is currently difficult as discussed in §1. 

To determine the bivariate PMF, tpp2, note that the 
number of protostars born in a time interval dt that will 
have final masses in the range dm f is 

d^Mp ^ jif*-)p{mf)dt dlnmf . (8) 
If we let tf{mf) be the time required to form a star of 
mass TO/, then integration of this equation gives the total 
number of protostars as 



••/ 



dlnmf ip{mf)tf{mf), 



= A{tf), 



(9) 
(10) 
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Fig. 1. The domain of integration for the PMF (cqs. 4, 5), 
At each value of the final protostellar mass, nif, integration over 
m extends from m = to m = m f . For eaeh value of the eurrent 
mass of the protostar, m, the integration over nif extends from 
rrif i = rnax(m, nif ) to the maximum mass of stars in the cluster, 
iiiu- In the plot, wc a.dopt irif = 0.033 and niu = 3. 

where we denote the average of some quantity x over the 
IMF as (a;) . Now, the characteristic accretion time scale 
for a protostar of mass m and final mass to/ is 



tacc{m,mf) 



dt 



din TO 



dm/dt 

so that equation (8) becomes 

d^Afp = AC'0(TO/)iacc dlnm dlnmf. 



(11) 



(12) 

Equations (3), (8) and (12) then give the final expression 
for the bivariate PMF, 



ipp2{m,mf) = 



■>P{mf%cc{m,mf) 

(tf) 



(13) 



This expression can be readily generalized to the case 
in which the accretion rate depends on more than two 
variables. 
The PMF (eq. 7) is then 



ipp{m) 



_ j_ 



ip{mf)tscc{'m,mf) dlnmf , (14) 



where mfj is defined in equation (6). The PMF is thus 
the IMF weighted by the accretion time, tacc = rn/rh, 
for all stars with final masses exceeding m. Note that 
as m. — > m.u the value of the integral approaches zero: 
there are very few protostars with masses close to the 
maximum value. 

3. ACCRETION HISTORIES 
3.1. Power-Law Accretion 

The PMF depends on both the accretion time, tacci 
and on the mean formation time, {tf) (eq. 14), so de- 
termining it requires evaluating the accretion histories of 
the protostars in the cluster. Several different models 
for protostellar accretion have been proposed, and in- 
deed measurement of the PMF would provide a powerful 
method for distinguishing among them. 

The most commonly used model for low-mass star for- 
mation is the inside-out collapse of a singular isothermal 
sphere (Shu 1977); we term this the Isothermal Sphere 



4 



McKee & Offner 



TABLE 1 
Symbol Key 



Symbol Definition 

ij^C Chabrier (2005) initial mass function 

m Instantaneous protostellar mass of a given star {Mq) 

ruf Final protostellar mass = initial stellar mass of a particular star (Mq) 

Lowest observable stellar mass (Mq) 

rriu Highest mass star that can be formed in the cluster (Mq) 

Afp Number of protostars 

A/* Star formation rate (number per unit time) 

tf Time to form a star of mass mf (yr) 

tfi Time to form a star of 1 solar mass without tapering (yr) 

m Instantaneous protostellar accretion rate (A'/0 yr^-') 

mi Final accretion rate of a 1 solar mass star without tapering (Mq yr""*^) 

ipp2 Bivariate protostellar mass function in terms of m and m f 

ipp Protostellar mass function in terms of m 



model. In this model, gas accretes from an isothermal 
gas in hydrostatic equilibrium (a "core") onto the pro- 
tostar at a constant rate determined by the temperature 
of the medium, 



mis = 1-54 X IQ-^Tj 



-6^7^3/2 



Mq yr- 



(15) 



where Ti = T/(10 K). This expression is based on the 
assumption that the gas is initially static. Hunter (1977) 
generalized the Shu solution to times prior to the initial 
formation of the protostar, so that at the time of proto- 
star formation the gas is in subsonic collapse. The most 
rapid collapse, at about 1/3 the sound speed, has an 
accretion rate 2.6 times greater than the Shu value. Fur- 
thermore, Li & Shu (1997) suggest that magnetic fields 
increase the accretion rate for low-mass protostars, typ- 
ically by about a factor 2. On the other hand, proto- 
stellar outflows are likely to eject some of the core mate- 
rial, reducing the final protostellar mass to a fraction 
ecore^^core, whcrc ecoro IS the corc star formation effi- 
ciency. Matzner & McKee (2000) estimated theoretically 
that ecore ~ 0.25 — 0.75; in a recent detailed study of 
cores in the Pipe Nebula, Rathborne et al. (2009) infer 
Ccorc = 0.22 ± 0.08 there. The effects of initial infall and 
magnetic fields on the one hand and protostellar out- 
flows on the other tend to cancel, but they render the 
estimate of the accretion rate in equation (15) somewhat 
uncertain. 

Since the time required to form a star via isothermal 
accretion scales directly as the mass of the star, this 
model does not work for high-mass star formation (Shu, 
Adams, & Lizano 1987). In order to treat high-mass 
star formation, McKee & Tan (2002, 2003) developed the 
Turbulent Core model as a generalization of the Isother- 
mal Sphere model. They considered a gravitationally 
bound clump of gas in which a cluster of stars is form- 
ing. Within this clump, individual stars form from bound 
cores. Both the clump and the embedded cores were as- 
sumed to be in approximate virial equilibrium, supported 
by internal turbulent motions. In terms of the surface 
density of the clump, Sci = Mc\/t^R% they found that 
the typical protostellar accretion rate is 



m ■ 



mxc 



mf 



3/4 



(16) 



where the coefficient mxc is proportional to the 3/8 
power of the pressure in the star-forming core and where 



the exponent j is related to the density profile of the core 
[p oc r~'^p) by 

(17) 

^ 2(3 -fcp)- 

To determine totc i they defined (j)p^ core as the ratio of 
the typical pressure of a star-forming core to the mean 
pressure of the clump. Since the pressure in a self- 
gravitating clump varies as S^j, their result corresponds 



to 

totc oc (j)p (,Qi.g5j(,j ■ They focused on high-mass star 
formation and allowed for the observed mass segregation 
of such stars by assuming that these stars formed on av- 
erage at about 30% of the half-mass radius of the clump; 
for this case, they estimated (pp. cote — 2. In the absence 
of such mass segregation, core — 1, and we adopt that 
value here. This reduces the accretion rate from their 
value by a modest amount [(1/2)^/^ = 0.77]. With this 
correction, the coefficient totc becomes 

mxc = 2.8 X 10-^J:1(^ Mq yr"!. (18) 

McKee & Tan (2003) set the remaining parameter in 
the accretion rate, j, by fixing the density power law at 
kp = ^ based on observations of clumps in which high- 
mass stars are forming. The precise value of kp is not 
important, however, since one can show that the value of 
mxc is within 10% of the value quoted in equation (18) 
over the range 1.3 < kp < 2. 

An assumption in the Turbulent Core model is that 
the turbulence is supersonic. McKee & Tan (2003) gave 
an approximate generalization of the accretion rate that 
includes the case in which the turbulence is subsonic and 
the accretion approaches the Isothermal Sphere value: 



^3/4 



m = mis 



2i 



1 + ^^ — ""f 



3/2 



where 



m/ / 



mxc 
rhis ' 



1/2 



(19) 



(20) 



For j ~ J, corresponding to fcp = 1, this is similar to the 
TNT model of Myers & Fuller (1992). 

A third model of accretion is Competitive Accretion, in 
which a group of stars in a common gravitational poten- 
tial accrete gas until it is exhausted or ejected (Bonnell 
et al. 1997, 2001a,b; Clark et al. 2007, 2009). These au- 
thors emphasize that the accretion rate depends on the 
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location of a protostar within the clump of gas and on 
the time evolution of the density in the clump. We take 
these effects into account indirectly by having the ac- 
cretion rate depend upon the final stellar mass so as to 
produce the prescribed IMF. Initially the accretion rate 
is governed by tidal effects, so that m oc m?/^; some 
of the stars fall to the center where they are virialized 
and then accrete at the Bondi-Hoyle rate (Bonnell et al. 
2001a), although the central, most massive star contin- 
ues to accrete at the tidal rate (Clark et al. 2009). We 
cannot take this change in the accretion rate into account 
in our model, Init we note that once the protostars are 
virialized, the mass accreted per dynamical time is small 
(Bonnell et al. 2001a; Krumholz et al. 2005). An im- 
portant feature of Competitive Accretion is that all the 
protostars cease accreting at the about the same time, 
presumably due to stellar feedback, as explicitly stated 
by Bonnell et al. (2001b) in their discussion of the IMF 
produced by Competitive Accretion. (Of course, in real- 
ity the gas removal is not instantaneous and the accretion 
will turn off more gradually; this can be treated by taper- 
ing the accretion rate, as discussed below.) Competitive 
Accretion is thus a constant time accretion model; this 
is in contrast to the Isothermal Sphere model, a con- 
stant rate accretion model in which the time for a star to 
form scales linearly with its mass. In order for all stars 
to form in the same time, the accretion rate must satisfy 
■m (xruf. The resulting model for Competitive Accretion 
has an accretion rate 



m \ 

TOCA =rni\ 



2/3 



m/. 



(21) 



where mi is the final accretion rate for star of unit mass. 
Bonnell et al. (2001a) show that the star-formation time, 
t/, is about equal to the initial free-fall time of the natal 
cloud, 

Myr, (22) 



tf c^ts= 0.435 fiH^i^ 



where nn, 4 = nH/(10^ cm~^) and fin is the mean density 
of hydrogen nuclei in the cloud. Integration of equation 
(21) shows that the characteristic accretion rate is related 
to tf by rhi = 3/t/. The principal approximation in our 
treatment of Competitive Accretion is the assumption 
that the star formation is in a steady state or is acceler- 
ating (§5). Our model thus applies to a cluster consisting 
of a number of small sub-clusters, each of which forms 
competitively, or to a sample of clusters. 

The Isothermal Sphere, Turbulent Core and Competi- 
tive Accretion models all fit the form 



TO = TOi I j nif 



3f 



(23) 



(see Table 2). For j < 1, this can be integrated to give 

m^-^ = (1 - j)mim/f-^ t. (24) 
The time to form the star (i.e., for to to reach to-j) is 

t/=t/iTO/i-^>, (25) 

where 



1 



(1 - j)mi 



(26) 



is the time to form a 1 Mq star. Note that if toi is 
expressed in units of Mq yr~^, for example, then tfi is in 



units of yr. The time scales for these models are included 
in Table 2, both the time to form a IMq star, i/(l), 
and the IMF-averaged formation time, {tf). (These time 
scales are for untapered accretion — see §3.2 below.) For 
the two-component Turbulent Core model (eq. 19), the 
star-formation time is 



tf 



(l + 7e^m/3/2)i/2 + i 



TO£ 
TOis' 



(27) 



which approaches the Isothermal Sphere value at low 
masses and the Turbulent Core value at high masses. 
The accretion time is then 



TO .N / ™ 

tacc = — = (1 - j) 

TO V"^/ 



1-i 



tf- 



(28) 



For protostellar masses significantly below the final stel- 
lar mass, the accretion time is longest for the Competi- 
tive Accretion model {j = |) and shortest for the Isother- 
mal Sphere model [j ~ 0). 

The models we consider here arc by no means exhaus- 
tive. It is also possible to contruct entirely different mod- 
els for the accretion history, for example;, those in which 
j > 1, so that the protostellar masses diverge at a fi- 
nite time (e.g., Behrend & Maeder 2001); such models 
require one to specify an initial protostellar mass, which 
is not well-defined. Myers (2009) has proposed a model 
that synthesizes the Isothermal Sphere and Competitive 
Accretion models. In this model, the protostar initially 
experiences freefall collapse where the accretion rate is 
comparable to rhg, while at late times the accretion rate 

approaches to* oc to*^^, which is close to the to* depen- 
dence of Bondi-Hoyle accretion. (Note that in our model- 
ing of Competitive Accretion, we have adopted the ml^^ 
dependence that Bonnell et al. 2001a found in their nu- 
merical simulations.) Models such as these that lead to 
an explosive growth in the stellar mass require an abrupt 
termination of the accretion. 

3.2. Tapered Accretion 

The accretion models presented above share one un- 
physical characteristic: the accretion drops to zero dis- 
continuously when the protostar reaches its final mass. 
Models such as these that lead to an explosive growth 
in the stellar mass require an abrupt termination of the 
accretion. Myers et al. (1998) addressed this problem 
by assuming that the accretion rate declined exponen- 
tially with time, to oc exp — (t/td)- A difficulty with this 
model is that the accretion continues indefinitely. More 
recently, Myers (2008) has included a model for a time- 
dependent dispersal of the protostellar core by protostel- 
lar outfiows. If the dispersal time is sufficiently short, 
the protostellar mass converges to a well-defined value. 

Even in the absence of prototellar outflows, one expects 
the accretion rate to taper off continuously. Observation- 
ally, Fedele et al. (2009) find that accretion falls below 
lO~"'^"'^M0yr~-'^ for nearly all stars by 10 Myr. However, 
the stars gain a negligible amoimt of mass via accretion 
after only a couple of Myr. For expansion wave solutions 
of the type in the Shu solution and the Turbulent Core 
model, the self-similarity of the accretion is broken when 
the expansion wave reaches the final mass, to/. If the 
core is immersed in a medium of uniform pressure, the ex- 
pansion wave will be reflected as a compression wave, and 
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TABLE 2 

Accretion Models 



Model 



Jf 



mi (Mq yr-i) tfi (Myr) (tf) (Myr) 



Isothermal Sphere 

Turbulent Core 1/2 3/4 

Competitive Accretion 2/3 1 



1.54 X lO-^Tf 



4.9 X lO-'^E^^ _^ 
6.9 X lO-^nJj'^4 



0.65 
0.40 S" 



-3/2 
-3/4 



0.25 Tj; 
0.29 E" 



-3/2 

-3/4 



= 0.435n, 



-1/2 



(*/> = % 



"Here _i = E(.i/{0.1 g cm ^) and nH,4 
scales are for untapered accretion. 

the accretion rate will decrease when that wave reaches 
the protostar (Stahler, Shu, & Taam 1980; McLaughlin 
& Pudritz 1997). 

A simple way to incorporate the decrease in accretion 
due to dispersal and boundary elfects is to reduce the 
accretion rate by a tapering factor [1 — (i/ij)"]: 



: nH/(10* cm Accretion rates and time 











1 - 









(29) 



where n > and where rhi is the final accretion rate for 
a 1 Mq star in the absence of tapering. Integration of 
this relation gives 



(1 — j)rhimf 



1 - 



1 



1 



The protostar reaches its final mass at a time 



n - 



1 



n 

n + l 
n 



V 



-J/ 



(1 - j)TOl ' 

t/ (untapered), 



t. (30) 

(31) 
(32) 



where it must be kept in mind that n > 0. 

We shall focus on the case n = 1, for which the forma- 
tion time is doubled over the untapered value given in 
§3.1. We note that the exponential tapering factor used 
by Myers et al. (1998) can be approximated by the n—1 
case for early times. In this case, one can solve for t{m), 
obtaining 

tim)= ; \ ' (n-1). (33) 

1 + [1 — (m/TO/)-^ JJJ-/^ 

Evaluation of 1 — t/t/ allows one to express the accretion 
rate in terms of only the masses. Defining 



untapered, n — 0, 

1 tapered, n = 1, 



(34) 



we can express the accretion rate for both the tapered 
and untapered cases as 

TOl I I m/f 1 - Snl 



(35) 
(36) 



The star-formation time for these two cases is 
tf = {I + (5„i)t/ (untapered). 

3.3. Mean Protostellar Mass 

Having described the accretion histories of the proto- 
stars, it is possible to infer the mean protostellar mass. 



(to) 



rmpp{'m)dh\ to. 



(37) 




Fig. 2. — Average values of powers of the protostellar mass, 
(m^), weighted by the Chabrier (2005) IMF for different powers of 
X clS cl function of the upper mass limit in the cluster, niu- 

Note that this is distinct from (to/), the final proto- 
stellar mass averaged over the IMF. However, because 
the PMF is given in terms of an integral (eq. 14), it is 
more straightforward to evaluate (to)^ using the bivari- 
ate PMF, 



din TO/ / dlnm 'mil)p2{m,mf) (38) 
Jo 



dlnm 



TO' 



3-Jf 



f- 



mi{tf) 



din J 



1 - Snl (to/to/) 



1-J 



1/2 



(39) 



The average star-formation time that enters this ex- 
pression is 

{tf) = {mf'-^Oi^ + Sni)tfi, (40) 

from equations (25) and (36). For a non-tapered accre- 
tion history {5ni = 0), the expression for the average 
protostellar mass reduces to 



I- J 



V 



(41) 



-jj {mf^ j) 

which must be evaluated numerically. In Figure 2, we 
have plotted (to=^) for the Chabrier (2005) IMF (see §4.2) 
for different exponents as a function of maximum mass in 
the cluster to facilitate the evaluation of equation (41). 

4. RESULTS FOR THE STEADY-STATE PMF 
4.1. Analytic Results for a Salpeter IMF 
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In order to see how the different accretion histories af- 
fect the PMF, we first consider a Salpeter IMF (Salpeter 
1955), 



1.35 



V'Sai(m/) ~ 1.35 — {mi<mf < rUy,), (42) 

where we have neglected the factor {rai/mu)^''^^ com- 
pared to unity in the normalization. For the case 
in which m„ is very large, the average mass (m/) = 
(1.35/0.35)m^; thus, to get an average mass of O.SM©, 
for example, requires = O.13M0. To compare with 
results from the truncated Chabrier IMF below, we also 
consider the case in which niu = ZMq and (m/) = 
O.4M0; this requires = O.I6M0. 

Since we are looking for simple analytic results in this 
subsection, we focus on the case of untapered accretion. 
The average star-formation time is then (eq. 25) 



where 



{ti) = 



1.35m/-J/0(j/) 
0.35 + j/ 



me 



(43) 



(44) 



For Isothermal Sphere accretion, (tf) is just {mf)tfi, or 
about half the time to form a one solar mass star if the 
mean mass is about O.5M0. 

With the aid of equation (23), the PMF in equation 
(14) becomes 



m /""" Vsai(?n/) 



1.35m/-^''^mi-J /""" din 



where 

and 



(1 - i)(0.35 + jf)mi°-^^+^fm^ 
a^{jf)mf^e°' 

a = 1.35+ jf -j 

/(m) ^ 1 



m 

TO,, 



(45) 
(46) 

- /(to/,,), (47) 
(48) 
(49) 



For masses large enough to be included in the IMF (to > 
m,), we have to/,, = m, so that 



tpp{m) oc m" 



(50) 



There are several points to note about this result for 
the PMF. First, since it is weighted by the accretion time, 
it is much flatter in the mass range to, < to to„ for 
Isothermal Sphere accretion (</> oc to~°'^^) than for Tur- 
bulent Core accretion {ip oc to~^'^) or Competitive Ac- 
cretion {ip oc to~^'^^); indeed, since all stars have the 
same formation time in the latter model, it just mimics 
the original Salpeter IMF. The Turbulent Core and Com- 
petitive Accretion models are shifted to lower masses 
compared to the Isothermal Sphere accretion model since 
they have longer accretion times at low mass, as shown 
in equation (28). 

Second, the PMF becomes depleted as to — > to„, since 
/(to) there. This occurs because protostars with 



final masses close to m„ spend most of their lives at lower 
masses as they grow by accretion; only a small number 
of protostars are actually in the final stages of growth. 

Third, the PMF is independent of the overall rate of 
star formation. This was clear from our general expres- 
sion for the PMF (eq. 14), in which tpp oc tacc/{tf) de- 
pends on the ratio of two star-formation times. 

The final point to note is that the coefficient in 'ipp{m) 
can be small, particularly as j 1; thus, most of the 
protostars can be at low masses. Evaluation of the mean 
protostellar mass using equation (41) shows that it is 
much less for the Turbulent Core and Competitive Accre- 
tion models than for the Isothermal Sphere model. For 
example, consider the case in which to, = O.I6M0 and 
TO„ = 3M0; we find that the average protostellar mass 
is (0.38, 0.16, 0.10) Mq for Isothermal Sphere accretion, 
the Tmbulcnt Core model and Competitive Accretion, 
respectively. 

4.2. Results for the Chabrier IMF 

The Salpeter IMF has the benefit of being easily inte- 
grable, but it is inaccurate at low masses. In this section 
we derive the mean mass for each accretion model using 
the Chabrier (2005) IMF: 



V'c = V'l exp - 



V'C = V'2TO 



(log TO — log 0.2) 
2 X 0.552 

(m; < TO < IM0 
-^•^^ (to > IM0). 



(51) 
(52) 



The Chabrier IMF is log-normal below 1 Mq and 
Salpeter above. The constants are determined by con- 
tinuity and by enforcing: /™" ■0(to) d In to = 1. In this 

section, we adopt fiducial values of to, ~ 0.033 Mq^ and 
TO„ = 3.0 Mq, which yield tpi ~ 0.35 and tjj2 — 0.16 = 
V'c(l)- With these coefficients, tpc gives an average stel- 
lar mass (including brown dwarfs) of (to*) ~ 0.4 Mq 

The PMF derived in equation (7) can be expressed 
in terms of the Chabrier IMF, j, jf, and to„. The Two- 
Component Turbulent Core model also depends upon the 
ratio TZrh = totc/tos. For the numerical examples pre- 
sented in this paper, we set the star-formation times for 
the Isothermal Sphere and Turbulent Core models to be 
equal, which the data in Table 2 show corresponds to 
TZm = 3.6. TZrh exceeds unity because it is evaluated for 
IM0 stars; for this value of TZm, the ratio of the accre- 
tion rates tor the two models is close to unity for pro- 
tostars of typical mass. Estimates of the star-formation 
time have typically been based on the Class lifetime, 
where most of the mass accretion was assumed to take 
place. However, observations by Enoch et al. (2008) indi- 
cate that the protostellar luminosities in the Class and 
Class I phases arc not substantially different, suggesting 
that significant accretion continues through much of the 
Class I phase. For untapered accretion, the PMF (eq. 7) 

^ This limit is derived assuming a 0.1 Mq minimum observable 
envelope mass, which converts its gas into a star with an efficiency 
factor of e 1/3. Reducing mi to the theoretically expected value 
of about O.OO4M0 (Whitworth ct al. 2007) has a small effect on 
the mean mass (at most a factor 1.1 for the Competitive Accretion 
model) and the median mass (a factor 1.25 for the same model); 
the ratio of the mean to median thus drops by no more than a 
factor 1.15. 
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Fig. 3. — The PMF for the four models with untapered accretion 
where m„ = 3 A/q. The Chabrier IMF is given by the solid line. 
For the Two-Component Turbulent Core model, we adopt TZm = 
3.6. 

for the IS, TC and CA models is 

. /™" '4'c{'nif)mf^~^f dlnmf 

(53) 

where we used equations (23) and (25) to evahiate the the 
acceleration time (tacc 1/m) and the mean formation 
time that enter the PMF. For tapered accretion, a factor 
[1 — (m/mfY^^]^/'^ must be included in the integrand in 
the numerator (see eq. 35). We plot the PMF for the 
four models in Figure 3 assuming that to„ = 3Mq. The 
Isothermal Sphere PMF peaks significantly to the right of 
the other models. As a result, it predicts a higher fraction 
of relatively massive protostars in the distribution than 
any of the other models, including the Chabrier IMF. 
This is a direct consequence of the weighting by the star- 
formation time, in which more massive stars have the 
longest accretion times in the Isothermal Sphere case. 
The Competitive Accretion and Turbulent Core PMFs 
contain fewer massive protostars than the Isothermal 
Sphere PMF because protostars in those models accrete 
more rapidly as they approach their final mass, thereby 
reducing the time they spend in the PMF. 

The mean protostellar mass, {m)^, is given by equation 
(41). With the Chabrier IMF, the Isothermal Sphere, 
Turbulent Core, and Competitive Accretion models give 
(m)p = 0.47, 0.16, O.O9M0, respectively, assuming that 
m„ = 3 Mq. The Isothermal Sphere mean value is 
larger than that derived in the previous section using 
the Salpeter IMF with mi — O.I6M0 since the Chabrier 
IMF turns over at O.2M0, so that the lower masses con- 
tribute less weight to the mean than in the case of the 
Salpeter IMF. The ordering of the values of the mean 
mass for these three models follows from the ordering of 
the values of the accretion time (eq. 28), since iacc is 
the weighting factor that enters the PMF (eq. 14) and it 
is largest at small masses for the Competitive Accretion 
model and smallest for such masses for the Isothermal 
Sphere model. 

Figure 5 shows {m)^^ as a function of the maximum 
cluster mass. The Isothermal Sphere model has a sig- 
nificantly higher mean protostellar mass than the stellar 



Fig. 4. — The PMF for the four models with untapered accre- 
tion(dotted), tapered accretion (dashed), and untapered accretion 
with accelerating star formation (dot-dashed), where rriu = 3 Mq. 
The Chabrier IMF is given by the solid line. For the Two- 
Component Turbulent Core model, we adopt 7?.m = 3.6. For the 
accelerating star formation models we use t = 1 Myr. 
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Fig. 5. — The mean mass (eq. 41) of the PMF, (m) ^, for the 
four models as a function of rriu in the untapered, non-accelerating 
case. The Chabrier IMF mean is shown by the solid line. For the 
Two-Component Turbulent Core model we adopt = 3.6. We 
restrict the Isothermal Sphere curve to < 5 Mq. 

IMF as a result of the long formation times of more mas- 
sive stars. In the Competitive Accretion and Turbulent 
Core models, the accretion rate accelerates with time so 
that protostars spend a larger fraction of their lifetime at 
low masses than in the Isothermal Sphere case, in which 
the accretion rate is independent of mass. Note that 
we assume that the Isothermal Sphere model is suitable 
only for stars below 5 Mq. When the accretion rate 
is modulated by turbulence, as in the Two-Component 
Turbulent Core model, the mean mass rises less steeply 
than for a pure isothermal sphere as a function of cluster 
mass upper limit. 

Figure 6 shows the ratio of the median to the mean 
protostellar mass as a function of the maximum stellar 
mass in the cluster, r7i„, for each of the accretion mod- 
els. This ratio decreases with to„ since the mean mass 
increases with m„ whereas the median mass is relatively 
insensitive to it. Since larger clusters can sample the 
rarer, more massive stars in the IMF, we expect large 
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Fig. 6. — The ratio of the median to mean protostellar mass 
of the PMF for the four models as a function of with non- 
accelerating untapered accretion. We adopt TZ^ = 3.6 model. The 
value for the Chabrier IMF is indicated by the solid line. Since 
the Isothermal Sphere model is primarily for low-mass stars, we 
restrict it to m,u < 5 Mq. 

clusters to have small ratios of the median to mean pro- 
tostellar masses, particularly for the Competitive Accre- 
tion and Turbulent Core models. Furthermore, since the 
accretion luminosity is proportional to the mass and the 
accretion rate, this graph suggests that the Protostellar 
Luminosity Functions will be quite different for the dif- 
ferent models, although we expect the differences to be 
less for tapered accretion (Paper II). 



The maximum expected protostellar mass, m 
a cluster of Ap protostars is given by 



p, max? 



1 



dlnmfTpp{mf) 



(54) 



The resulting values of nip^ max are portrayed in Figure 
7. The plot illustrates that the expected value of the 
maximum protostellar mass, TOp,max, for the Competi- 
tive Accretion and Turbulent Core models with A/"p ^ 100 
stars can be significantly less than the maximum stellar 
mass, rriu, which is determined by the Class II stars in 
the cluster that obey a normal IMF. 

5. ACCELERATING STAR FORMATION 

From an analysis of the pre-main-sequence stars in the 
Orion Nebula Cluster, PaUa & Stabler (1999) concluded 
that the star formation there has been accelerating. Palla 
& Stabler (2000) extended this analysis to seven other 
star-forming regions and found evidence for acceleration 
in all but one case. They attributed this acceleration 
to contraction of the parent molecular cloud, which they 
surmised was a quasi-static process. They inferred ex- 
ponentiation times ranging from 1.0 Myr for p Oph and 
IC 348 to 3.3 Myr for NGC 2264. The' shortest of these 
acceleration times is not that much greater than the typ- 
ical star-formation time of 0.54 Myr found by Evans et 
al. (2009). Here we determine the effect of accelerated 
star formation on the PMF. 

The evolution of the PMF in time is governed by a 
continuity equation. To conform with standard practice, 
we define n(m, mj, t)dmdmf as the number of protostars 
in the mass range dm and final mass range dm f at time 
t. Under the assumption that the stars are born with an 



Fig. 7. — The maximum protostellar mass, nip^ max as a function 
of the number of protostars, A/'p in the untapered, non-accelerating 
case. The maximum stellar mass is rriu = 3 Mq. We adopt TZm = 
3.6 for the Two-Component Turbulent Core model. 

initial mass function tfj{mf), the continuity equation for 
n is then 



at am m f 



(55) 



where AC (t) is the rate at which stars are born at time 
t, S{x) is the delta function, and the factor mf~^ al- 
lows from the conversion from n dm f to ^ dlnmf. The 
Green's function for this problem, which we denote by 
G{t - io), is the solution for J\f^{t) ^ d{t - to). Writ- 
ing the protostellar mass as an explicit function of time, 
m = /i(i — to), we have 

G{t - to) = S[m ~ ii(t - to)\ H{t - to), (56) 

mf 

where H{t — to) is the Heaviside step function. The gen- 
eral solution is then 



— oo 

t 



G{t ~to)M*,odto, 

S[m - ii{t ~ to)] ^^^^^ a; odto. 
mf 



(57) 



(58) 



Note that t — to is the age of a star at time t that was 
born at time to- Let t^ be the age of a star of mass m 
and final mass mf. We can then rewrite the (5- function 
as 

S[m-^,{t-to)] = ^-^^^^^^^, (59) 

m 

where m = dfi/dt is the accretion rate. The solution is 
then 

1p{mf)J\f^{t ~tjn) 

n{m,mf ,t) = . 

mfm 

To convert this to the bivariate PMF, note that 
n{m, mf, t) mmf 



(60) 



il)p2{m,mf,t) 
so that 



/X" ^"^/ /o"^ '^"^ '^("^' "^Z' ^) ' 



7/'(m/)— 7V;(t-t„) 
m 



!Zt dlnmfiPimf) /^^ dt^N*{t - tr 



(61) 



(62) 
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In a steady state [^(t — tm) = eonst], this reduees to 
the result given in equation (13). The protosteUar mass 
function for an arbitrary star formation history is then 
obtained by inserting this result into equation (7). 

As a simple model of accelerating star formation, we 
assume an exponentially increasing birthrate, 

M,(i^t™) =^;,oe(*-*'")/^ (63) 

where AC^o = M*{t = 0) is the current birthrate. Substi- 
tuting into equation (62), we find 

fp{mf)ti,cc exp -{tm/r) 



V'p2 = 

The PMF is then 

ii)p{m) = 



t(1 - exp{-tf/T)) 



(64) 



r(l - exp{-tf/T)) 

V' (m/ )tacc e In m/ , 



(65) 



which reduces to equation (7) for r — >■ oo. 

The PMF depends on the time scale t„i, the age of a 
protostar of mass m and final mass nif. If the accretion 
is untapered, then the protostellar mass grows according 
to equation (24) and tm is given by 



tm = tfl 



m 



(66) 



For the Two-Component Turbulent Core model, tm is 
given by 



_ 2 ^Ttjmfy^m + l-l 

The value of tm for the case of tapered accretion has been 
given in equation (33). Figure 4 shows ipp{m) in the ac- 
celerating case with untapered accretion. We consider 
the case with r = 1 Myr, which is approximately twice 
the average star-formation time. In comparison to the 
steady star formation (dotted lines), the PMF peaks are 
shifted towards lower masses for accelerating star forma- 
tion except in the Competitive Accretion case, where the 
peak mass increases from ~ O.O5M0 to ^ O.OSMq. 

6. COMPARISON WITH OBSERVATIONS 
6.1. The Star Formation Timescale 

As discussed in the Introduction, it is not presently 
possible to directly measure the masses of the proto- 
stars in a star-forming region. The protostellar lumi- 
nosity function is subject to direct observation and will 
be discussed in Paper II. Nonetheless, it is possible to 
carry out an approximate comparison between observa- 
tion and theory by comparing the average star-formation 
time observed in low-mass star-forming regions with the 
theoretical values predicted in §3. 

First, we show that the average observed star- 
formation time as determined by number counts is 
the IMF-averaged formation time, (tf), not the PMF- 
averaged formation time, {tf)^. Evans et al. (2009) deter- 
mined the average star-formation time, (i/)obs, by com- 
paring the number of protostars with the number of Class 
II sources: 



(^/)o 



(in) Ml' 



where Mi and (tn) arc the number and average lifetime 
of Class II sources. Although nothing is known at present 
about the mass dependence of the Class II lifetime, we 
allow for the possibility that there is such a dependence 
by using the IMF-averaged value, (tu), in this equation. 
The number of protostars is just the star-formation rate 
times the IMF-averaged lifetime (eq. 10), and similarly, 
the number of Class II sources is just Mi = M*(tii). As 
a result, we have 



(^/)obs ^ M(^/) ^ jtf) 

(ill) A{tu) {tu}' 



so that 



(*/)obs = (*/)• 



(69) 



(70) 



bs 



(68) 



Thus, the mean formation time inferred from the ratio 
of the number of protostars to the number of Class II 
sources is equal to the IMF-averaged value of the for- 
mation time. By contrast, if there were a method of 
determining the formation time of each observed proto- 
star, the average of these times, ~ J dliimijjp{m)tf , 
would be quite different. We generalize this expression to 
the case of accelerating star formation in the Appendix. 

Evans et al. (2009) report an average star-formation 
time, (tf), of 0.54 Myr after correcting for extinction for 
five local star-forming regions. This value is the sum 
of the estimated Class and Class I lifetimes and as- 
sumes a Class II lifetime of 2 Myr. Evans et al. (2009) 
find that the lifetimes vary significantly over their sam- 
ple of clouds, where the lowest lifetimes correspond to 
the smallest clouds containing the fewest protostars. The 
lifetimes have an uncertainty of order ±0.1 Myr. In our 
comparison, we focus on Perseus, Ophiuchus, and Ser- 
pens, which are much more massive than Lupus and Cha 
II. They have a combined average star-formation time of 
0.56 Myr. 

As discussed in the Introduction, by adopting the pro- 
tostellar lifetimes from Evans et al. (2009), we implic- 
itly assume a direct mapping between Class 0, I and 
Stage 0, I. The core mass estimates obtained for these 
regions by Enoch et al. (2009) support this assumption 
for both Perseus and Serpens since all but one Class I 
envelope exceed O.IMq in Perseus and none in Serpens. 
In contrast, Enoch et al. (2009) report that about half 
of the p Ophiuchus Class I sources have envelopes with 
masses < O.IMq, so that they are not true protostars 
by our definition. This suggests that the value of {tf)^^,^ 
we adopt from Evans et al. (2009), which includes the 
Class I sources with envelopes less than O.IMq, overesti- 
mates the mean protostellar lifetime in that region. Since 
Ophiuchus already has the shortest lifetime of the three 
clusters we consider, it is possible that the rate of star 
formation in this cluster has recently slowed down. Bear- 
ing in mind these caveats, in our comparison we adopt 
conservative error estimates that are comparable to the 
error resulting from the inclusion of the small-envelope 
sources. 

To compare with the observations, we use the more 
evolved Class II sources to estimate the cluster upper 
mass limit, m„. Since the Class II sources are about 
three times more numeroTis than the Class and Class 
I objects, their luminosities provide a reasonable expec- 
tation for the most massive star forming in the cluster. 
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TABLE 3 
Star-Forming Region Properties 



Region T (K)=' S (g cm-2)b fin (10"* cm-=^)'= (to,i> (Myr)<^ 



Perseus 10-13 0.06 1.0 0.72 

Serpens 10-15 0.06 0.9 0.56 

Ophiuchus 12-20 0.08 1.4 0.40 



Average 007 LI 0.56 

^Tlic citations for the tcnipcratures in each of tlie regions are Foster et 
al. (2009); Enoch et aL (2008); Andre et al. (2007) 

''Average gas column of the clumps with > 10, where the clumps 
are assumed to be spherical. Masses and sizes are supplied by Evans 
et al. (2007) and Enoch et al. (2007). The gas with Ay > 10 ap- 
proximately corresponds to the minimum column density for star for- 
mation: JV(H2) = 8 X 10^1 cm-2 (Onishi et al. 1998), where Ay = 
Af(H2)/(0.94 X 10^1 cm-2). 

"^The densities are calculated by identifying clumps of gas with Ay > 
10. The average density of each region is calculated assuming spherical 
symmetry. The mass weighted average of the clumps is reported here. 
Data are derived from Evans et al. (2007) and Enoch et al. (2007). 

''Average extinction corrected (Class + Class I) lifetimes (Evans et 
al. 2009) 



Evans ct al. (2009) report a maximum extinction cor- 
rected Class II luminosity of ~ 90Lq, suggesting a max- 
imum stellar mass of ~ 3M© (Palla & Stabler 1999). As 
a counterpoint, if we assume that the combined popula- 
tion of Class II sources is drawn from the Chabrier IMF 
with maximum stellar mass of 120 Mq then a statistical 
argument suggests that 12 of the Class II objects should 
have masses greater than 3 M©. Perseus, containing 244 
Class II sources, should have five such stars. This differ- 
ence between the observed IMF and the expected one is 
partially due to selection, since the clusters were chosen 
to have only low-mass stars. It is also possible that the 
conditions in these clouds work against the formation of 
massive stars. 

We use the model definitions from Table 2 to infer the 
physical parameters for each model given an average star- 
formation time of {tf) = 0.56 Myr. We give the observed 
physical parameters values for local molecular clouds in 
Table 3. Each of the models constrains a different physi- 
cal parameter (Table 4). For the Competitive Accretion 
model, the relevant star-formation time, (tf), is about 
equal to the freefall time of the entire clump, so that the 
lifetime depends only upon the average density. For the 
Two-Component Turbulent Core model, we assume that 
the thermal and turbulent contributions to the accretion 
are comparable by setting TZm = 3.6 (see §4.2). (The 
Turbulent Core model was developed for high-mass star 
formation, and the Turbulent Core contribution grows in 
importance as the stellar mass increases, eq. 19.) 

To compare the models with accelerating star forma- 
tion with observation, we first use a least-squares ap- 
proach to re-analyze the data assembled by Palla & 
Stabler (2000), who fit their data "by eye." We find 
T = 0.9 Myr for Ophiuchus, in good agreement with 
the Palla & Stabler (2000) value of r = 1 Myr. At face 
value, this suggests that star formation is rapidly accel- 
erating. However, the ages of stars included in the fit are 
greater than 0.5 Myr, so that one cannot use this trend to 
reliably infer information about more recent star forma- 
tion activity. To emphasize this point, observations by 
Enoch et al. (2009) find only 3 Class sources, which, 
when compared to the relatively larger numbers of Class 



I sources, suggest that current star formation in Ophi- 
uchus is likely deceleraiing. A fit of IC348, an active; rc;- 
gion of Perseus, gives r = 2.2 Myr rather than 1 Myr as 
reported by Palla & Stabler (2000) . Several other regions 
discussed by Palla & Stabler (2000) show a turnover in 
the number of stars at recent times, further highlighting 
the point that the star formation rate may not be mono- 
tonic and is not necessarily well represented by a single 
exponential. Despite this caveat, we select r = 2 Myr 
as a fiducial value. For Perseus and Ophiuchus, we indi- 
vidually adopt T = 2 Myr and r = 1 Myr, respectively. 
We adopt 2 Myr for the Class II lifetime, which is the 
same value Evans et al. (2009) use to estimate the ob- 
served protostellar lifetime. Note that it is possible that 
the Class II lifetime depends on stellar mass, but in the 
absence of any relevant data we neglect this possibility 
here. 

The results summarized in Tables 3-5 show the inferred 
physical conditions in an average star-forming cloud (like 
Serpens), in Perseus and in Ophiuchus for each model 
for the cases of tapered and untapered accretion, ac- 
celerating and non-accelerating star formation. The in- 
ferred temperatures for the Isothermal Sphere and Two- 
Component Turbulent Core models are closer to the ob- 
served values when tapering and/or acceleration is al- 
lowed for. Tapering or allowing for acceleration gives 
a good fit for both the Turbulent Core model and the 
Competitive Accretion model in Perseus. In Serpens, 
non-accelerating star formation with untapered or ta- 
pered accretion gives good agreement for these models. 
In Ophiuchus, both models do best for the untapered, 
non-accelerating case. 

However, as discussed in §3 above, the accretion rates 
for the various models are uncertain for several reaons: 
they do not allow for the effects of protostellar winds in 
disrupting the protostellar cores, they do not allow for an 
initial infall velocity, and they do not include the effects 
of magnetic fields. As a result of these competing effects, 
the fiducial accretion rate coefficients in Table 2 may be 
either higher or lower than those we have adopted, thus 
impacting our estimates of the physical parameters. If we 
assume that this uncertainty is a factor two in the accre- 
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tion rate, then this corresponds to a factor 2^/"^ 1.6 in 
the inferred temperature in the Isothermal Sphere model, 
a factor 2"*/^ ~ 2.5 in the inferred column density in the 
Turbulent Core model, and a factor 2^ — Am the inferred 
density in the Competitive Accretion model. It must also 
be borne in mind that there are uncertainties in the ob- 
served values of the physical parameters describing the 
clouds, particularly the lifetimes. With allowance for an 
overall factor of two uncertainty in the accretion rate rel- 
ative to the actual values, we summarize the consistency 
between the models and observed regions in Table 7. We 
find that the Isothermal Sphere model is consistent with 
the data for both Perseus and Ophiuchus for all but the 
non-accelerating, untapered model. One might question 
this consistency, given that the temperature is known to 
much better than the factor 1.6 we are allowing; how- 
ever, we are allowing for uncertainty in the coefficient 
of T'^/^ in equation (15), not for errors in the temper- 
ature. The Two-Component Turbulent Core model is 
consistent with the data for all but the non-accelerating, 
untapered model; it would have the same level of consis- 
tency with the data as the IS model if TZm were lowered 
from 3.6 to 1.7. (Based on eq. 27, this corresponds to 
increasing the value of mj for which i/ is closer to the 
TC case than to the IS case from about 0.5 Mq to about 
1.2 M0.) The Turbulent Core and Competitive Accre- 
tion models are consistent with the Perseus data for all 
but the non-accelerating, tapered model, but they are 
consistent with the Serpens and Ophiuchus data only for 
the non-accelerating models. 

We note that the free-fall time corresponding to the 
mean density in the regions we have considered is ~ 
0.4 Myr, which is comparable to the mean star formation 
times (which is consistent with Competitive Accretion), 
but short compared to the acceleration time. Concep- 
tually, the Competitive Accretion model involves rapid 
acceleration, but numerical simulations are needed to de- 
termine whether the observed relatively long acceleration 
times are consistent with the model. 

Since a number of the models predict reasonable phys- 
ical parameters, uncertainty in the mean temperatures, 
densities, columns, and lifetimes contributes to our in- 
ability to determine which models are consistent with 
observation. For comparison in the context of these un- 
certainties, we plot the physical parameters of each of the 
model as a function of (tobs) in Figure 8. This plot also 
illustrates that in all cases the inferred temperature, col- 
umn density, and density increase with decreasing star- 
formation times. 

6.2. Discussion 

In this section, we address various sources of uncer- 
tainty inherent in the calculations. As described in §3.1, 
the instantaneous accretion rate may be larger than es- 
timated at early times and may fluctuate strongly on 
short timescales. We have also neglected possible tem- 
poral variations in the core star formation efficiency-i.e., 
the fraction of the mass in the core that actually accretes 
onto the star-in modeling the protostellar accretion rate. 
Finally, the values of n, (tu) and r have some intrinsic 
uncertainty and may vary between clouds. We discuss 
these three parameters below. 

Tapering of the accretion is very plausible, since we 
expect accretion to decrease as the infall phase ends and 
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Fig. 8. — The physical parameter for each model versus the ob- 
served star-formation time, {tf)^^^^■ We assume n = 1, r = 1 Myr, 
and (til) = 2 Myr in the tapered and accelerating cases. Ophiuchus 
(left), Serpens (middle), and Perseus (right) are overlaid with hor- 
izontal error bars for the uncertainty in the measurements and 
vertical error bars for the uncertainty due to the model accretion 
rates. 

the core mass is depleted. Measurements of the proto- 
stellar luminosity appear inconsistent with constant or 
increasing accretion through the Class I phase, a topic 
we shall address in Paper II. However, there are no direct 
estimates of n, and the exact tapering function is poorly 
constrained by observation. Unfortunately, we find that 
we cannot definitively constrain n on the basis of the 
observed star-formation time due to variation between 
clouds and uncertainties in both the physical parameters 
and the accretion models, as discussed above. 

Evans et al. (2009) estimate an uncertainty of ±1 Myr 
for the Class II lifetime, a variation of 50%. As shown 
by equation (68), the observed star-formation time varies 
directly as (tu) in the non-accelerating case. The accre- 
tion rate then varies as l/(tn), so that T'^/^ cic l/(iii) in 
the Isothermal Sphere case, cx l/{tu) in the Tur- 

1/2 

bulent Core case, and oc l/(iii) in the Competitive 
Accretion case. Thus, a longer Class II lifetime would 
decrease the inferred temperature, column density and 
density accordingly. Our analysis in §6.1 suggests that 
a longer Class II phase would worsen agreement of the 
Isothermal Sphere and Two-Component Turbulent Core 
models with the observations, since the implied temper- 
atures are already somewhat low in most cases. Con- 
versely, a shorter Class II lifetime would improve agree- 
ment with observation for many of these models, with the 
exception of the tapered, accelerating Isothermal Sphere 
model for Ophiuchus. 

As discussed in §5, Palla & Stabler (2000) derived r 
values between 1 — 3 Myr for a collection of eight 
clouds. They performed the fits by eye so that their es- 
timates of the e-folding time are not precise. We used 
a least-squares fit to get improved values for the accel- 
eration times, but it must be borne in mind that the 
actual star formation histories may be much more com- 
plex than the simple exponential form that we have used 
to fit them. 
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TABLE 4 





Model Properties for (t/)oi,s 


= 0.56 Myr 








N o n- A c c 1 cr at i n g 


Accelerating''' 


Model 


i*aiaui(M cr 


L"iilai)('r(HL 


TaiKTc<r 


L 111 aix'rcd Tapcicd'* 


Isothermal Sphere 


T(K) 


5.9 


9.3 


7.8 12 


2C Turbulent Core 


T(K) 


3.9 


6.1 


5.4 8.6 


Turbulent Core (TC) 


S (g cm ^) 


0.041 


0.10 


0.083 0.21 


Competitive Accretion 


fiH (10"' cm~^) 


0.60 


2.4 


1.8 7.1 



''The fiducial (ty) is the same as {tf)^^^ for Serpens, so the data for Serpens should be 
compared with these values, 
br = 2 Myr; <t„) = 2 Myr 
"^n = 1 

''We fix TZ^ = 3.6. 



TABLE 5 

Model Properties for (tf)^^^^ = 0.72 Myr: Perseus 









Non- Accelerating 


Accelerating 


Model 


Parameter 


Observed Value 


Untapered Tapered*" 


Untapered Tapered*" 


Isothermal Sphere 


T (K) 


10-13 


5.0 7.9 


6.6 10 


2C Turbulent Core<^ 


T (K) 


10-13 


3.3 5.2 


4.6 7.3 


Turbulent Core (TC) 


S (g cm 


0.06 


0.029 0.074 


0.061 0.15 


Competitive Accretion 


riH (10"' cm-3) 


1.0 


0.37 1.5 


1.1 4.5 



*r = 2 Myr, the value for IC 348; (tu) = 2 Myr 
*>«= 1 

"^We fix TZ^ = 3.6. 



TABLE 6 

Model Properties for {t/>oj,s = O-^O Myr: Ophiuchus 









Non-Accelerating 


Accelerating ^ 


Model 


Parameter 


Observed Value 


Untapered Tapered*" 


Untapered Tapered*" 


Isothermal Sphere 


T (K) 


12-20 


7.4 12 


13 20 


2C Turbulent Core 


T(K) 


12-20 


4.8 7.7 


8.7 14 


Turbulent Core (TC) 


S (g cm-2) 


0.08 


0.064 0.16 


0.22 0.55 


Competitive Accretion 


riH (10* cm-3) 


1.4 


1.2 4.7 


7.4 30 



a-r = 1 Myr; (tu) = 2 Myr 
*"n = 1 

<^We fix TZrh = 3.6. 



TABLE 7 

Model St^biary 



Ophiuchus Serpens Perseus 

Model UN a TN UA TA UN TN UA TA UN TN UA TA 



Isothermal Sphere / / // // // o/ // 

2C Turbulent Core o / // o/ Z/ o/ // 

Turbulent Core (TC) / / 00//00///0 

Competitive Accretion / / 00//00///0 



''UN = Untapered Non-accelerating,TN = Tapered Non-accelerating, UA = Untapered Acceler- 
ating, TA = Tapered Accelerating 



Wc have analyzed the Protostellar Mass Function 
(PMF), which is the Present-Day Mass Function 
(PDMF) of a cluster of protostars. The PMF builds on 
the IMF, which measures the final mass distribution of a 
cluster of stars after completion of the formation process. 
The PMF depends on the mass dependence of the forma- 
tion time of the stars, tf, and therefore can in principle 
provide an observational diagnostic of theories of star 
formation. Direct observation of the PMF must await 



further observational and theoretical advances, since it 
is not currently possible to infer the mass for proto- 
stars because they are generally extremely obscured and 
the luminosity is dominated by accretion. Nonetheless, 
it is possible to compare theoretical models with mea- 
surements of the mean lifetime of protostars in different 
molecular clouds as we have done here. Furthermore, the 
PMF is the basis for the Protostellar Luminosity Func- 
tion, which is directly accessible to observation (Paper 
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Wc have made two key assumptions in inferring the 
PMF. First, since no stars above 3 Mq are seen in the 
low-mass clusters we have analyzed (Evans et al. 2009; 
Enoch et al. 2009), we have assumed that none of the 
smaller number of protostars in these clusters will grow 
to more than 3Mq; wc therefore imposed an upper cut- 
off of 3Mq on the IMF in our analysis. Given the num- 
ber of Class II sources, which have very nearly reached 
their final masses, we note that 12 Class II sources with 
m > 3Mq would be expected in these clusters based on a 
Chabrier (2005) IMF. These clusters were selected on the 
basis of their proximity so as to allow a thorough study 
of the Class 0, I, and II sources, so it is unlikely that 
this anomalous IMF is due to a selection effect. The sec- 
ond key assumption we have made is that the accretion 
rate is a simple function of only the current protostel- 
lar mass, the final mass and the time. Wc thus do not 
allow for variations in the accretion rate due to a brief 
high-accretion Larson-Penston phase or to temporal fluc- 
tuations in the accretion rate (although insofar as such 
fluctuations arc random and there is a statistically large 
sample of protostars, they should not significantly affect 
the PMF). 

We consider four accretion rate histories: the classi- 
cal Isothermal Sphere accretion (Shu 1977), the Turbu- 
lent Core model (McKec & Tan 2002, 2003), a blend of 
the two (Two-Component Turbulent Core, 2C Turbulent 
Core), and an analytic approximation for the Competi- 
tive Accretion model (Bonnell et al. 1997, 2001a). There 
are substantial uncertainties in the accretion rates for 
each model: In all cases, one must allow for the effect 
of protostcUar outflows, which can reduce the accretion 
rate by a factor of a few (Matzner & McKee 2000). For 
the first three, there is a countervailing correction needed 
to allow for an initial infall velocity. Our approximation 
for the Competitive Accretion model captures many of 
its essential features, but since the model itself is based 
primarily on numerical simulations, there is no fully an- 
alytic form for it. In comparing the models with obser- 
vation, we assume that the star formation is steady or 
accelerating; since the Competitive Accretion model has 
been developed for the evolution of individual star clus- 
ters, the comparison with observation is valid for this 
model only if a number of clusters are sampled, either 
because a forming cluster is comprised of a number of 
sub-clusters or because data from different clusters are 
averaged together. 

The mean protostcUar mass (Fig. 5) and the ratio of 
the median mass to the mean mass (Fig. 6) depend sensi- 
tively on the accretion history. The Turbulent Core and 
Competitive Accretion Models have accretion rates that 
increase with mass and therefore with time (to oc to-', 
with j = i, I for the two models respectively). As 
a result, protostars of a given final mass, m/, spend a 
smaller fraction of their lives at high mass than in the 
Isothermal Sphere model. Furthermore, these two mod- 
els have accretion rates that increase with mj, so that 
it takes less time to form a high-mass star than in the 
Isothermal Sphere model. Both effects are stronger for 
the Competive Accretion model. As a result, the mean 
protostcUar mass increases systematically from the Com- 
petitive Accretion model to the Turbulent Core model 



to the Two-Component Turbulent Core model to the 
Isothermal Sphere model. The ratio of the median to 
mean protostellar mass follows the same ordering, and 
the same effect shows up in the plots of the PMF in Fig- 
ure 3. 

A common feature of all the accretion models is that 
the accretion rate remains constant or (usually) increases 
until the time at which the protostar reaches its final 
mass, when it abruptly ceases. In reality, as pointed 
out by Myers et al. (1998), the accretion will turn off 
gradually. To allow for this, we have inserted a factor 
[1 — (t/t/)"] into the accretion rates; we refer to this as 
tapered accretion. In practice, we focused on the case 
n = 1 , which gives an accretion time t / twice as long as 
would be expected in the absence of tapering. This has 
the effect of increasing the temperature, column density 
or density of the model needed to match a given observed 
formation time. For example, in the Isothermal Sphere 
model, the accretion rate is proportional to T^/^, so the 
temperature needed to match the observations of a given 
formation time is 2^/"^ times greater for a tapered model 
than for an untapered one. Not only is tapering phys- 
ically plausible, it also generally results in models that 
are in better agreement with observation. As shown in 
Figure 4, tapering moves the peak of the PMF to higher 
masses since stars spend a larger fraction of their lives at 
high mass when the accretion slows down at the end of 
the accretion process. 

The rate of star formation should accelerate in time 
in a contracting gas cloud, and Palla & Stabler (2000) 
found direct evidence for such acceleration in a num- 
ber of nearby star-forming clusters. We generalized our 
analysis of the PMF to the case in which the star for- 
mation rate is time dependent in §5. For simplicity, we 
have assumed that the acceleration applied only to the 
rate at which stars formed, not to the accretion rates of 
individual stars. This is a reasonably good approxima- 
tion for the cases we analyzed, which have star- formation 
times that are significantly smaller than the time scale 
for acceleration. Moreover, the approximation is even 
better for for the Isothermal Sphere model, since the ac- 
cretion rate depends on the temperature and radiative 
losses maintain an approximately constant temperature. 
On the other hand, the time scale for acceleration is of- 
ten comparable to or less than the mean lifetime of Class 
II sources, so the ratio of the number of protostars to 
the number of Class II sources is larger than in the non- 
accelerating case. As a result, as shown in the Appendix, 
the "observed" star-formation time, which is given by 
equation (68), exceeds the actual star- formation time. 
We find that acceleration does not have a substantial ef- 
fect on the PMF. Rather, its primary effect is to reduce 
the inferred time scale for the formation of individual 
stars, thereby increasing the inferred temperature, col- 
umn density or mean density, depending on the accretion 
model. 

In the absence of any direct information on proto- 
stellar masses, we were able to carry out only a very 
crude comparison with observation: Using the observed 
star-formation time scales in two different clusters, we 
computed the implied temperature (Isothermal Sphere 
model), surface density (Turbulent Core model), and 
mean density (Competitive Accretion model), and then 
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compared with the observed values of these parameters. 
We found that the tapered accretion and accelerating 
star formation models were somewhat better than unta- 
pered, non-acceleration models, but we could not draw 
any firm conclusions due to uncertainties in both the ob- 
servations and in the models, which have accretion rates 
that arc probably uncertain by a factor of 2. In addition, 
the molecular clouds have an unknown internal structure 
and the IMF can have significant statistical and perhaps 
physical fluctuations from one cloud to another. In Pa- 
per II we shall show that the Protostellar Luminosity 



Function is a more powerful diagnostic for inferring the 
accretion mechanism. 
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APPENDIX 

OBSERVED STAR FORMATION TIME FOR ACCELERATING STAR FORMATION 
For a star formation rate that varies exponentially in time, the number of protostars is given by 



AAp(t = 0): 



dmn{m,mf,t = 0) 



-- / dnif / 

J mg_ J 



=7V;,or(l - e" 



(Al) 

(A2) 
(A3) 



In the time-dependent case, the number of Class II sources in the mass range druf is the number formed in the time 
interval tf < to < tf + tu, 

rtf+tu 

dMu{mf) = ijj{mf)dlnmf / Af*{t-to)dt. (A4) 

Jtf 

For an exponentially increasing star formation rate (eq. 63), the total number of Class II sources is then 



(A5) 



where the average is over the IMF. With the aid of equation (A3), equation (68) then implies that the mean observed 
star-formation time is 

_ (fa) (1 - e-*//-) 
W/obs - ^e-*//^ (1 - e-*"/^)) ■ ^ ' 

In the limit of steady star formation (r — > oo), this approaches the actual value of the average star-formation time, 
(tf). For tf <C r, which is true for most of the models we have considered, the observed star-formation time is related 
to the actual value by 

(tu) 



(tf) 



obs 



T (l - e-<*">/-^) 



{tf), 



(A7) 



where we have made the approximation (1 — exp{—tu/T)) ~ 1 — exp(— (iii)/r). The observed formation time, (i/obs), 
thus exceeds the actual value, {tf), since the number of Class II sources is suppressed. In the text we evaluate equation 
(A6) for tu = 2 Myr for each of the accretion cases considered. 
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